Generalized Numerical Renormalization Group for Dynamical Quantities 
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In this paper we introduce a new approach for calculating dynamical properties within the numerical 
renormalization group. It is demonstrated that the method previously used fails for the Anderson 
impurity in a magnetic field due to the absence of energy scale separation. The problem is solved 
by evaluating the Green function with respect to the reduced density matrix of the full system, 
leading to accurate spectra in agreement with the static magnetization. The new procedure (de- 
noted as DM-NRG) provides a unifying framework for calculating dynamics at any temperature and 
represents the correct extension of Wilson's original thermodynamic calculation. 



Quantum impurity models and their low-temperature 
properties are of central importance in condensed matter 
physics. They show characteristic many-body effects like 
the screening of a local moment by conduction electrons 
(the Kondo effect) which was first observed in measure- 
ments on dilute magnetic impurities in metals (see [0). 
More recently, artificial nanostructures (quantum dots 
1^ or surface atoms probed by STM j|,^) with tunable 
parameters provided new representations of the Ander- 
son or Kondo model (^J^. In theory, a very fruitful line 
of research was opened by the development of dynami- 
cal mean-field theory (DMFT) where correlated lat- 
tice systems are mapped onto effective impurity models 
which are then accessible in a controlled way |^ . 

In all the above areas, progress depends sensitively 
on the existence of a reliable calculational method that 
can provide static and dynamic (spectral) properties in 
the full energy range. Wilson's numerical renormaliza- 
tion group 1^ gave the first quantitative description of 
the Kondo effect. In systems with very different energy 
scales (small Kondo temperature Tk, large bandwidth) 
it is the only technique that can do so. In the original 
calculation Wilson focused on obtaining thermodynamic 
expectation values like the impurity susceptibility by it- 
erative diagonalization. Each iteration step was shown to 
correspond to a certain temperature where expectation 
values could be obtained with great precision. Later, 
the method was extended to zero temperature dynam- 
ical properties by several groups and applied to a vari- 
ety of proble ms [10|jTl| ], including recent DMFT calcula- 
tions (e.g. 13|). In these calculations the additional 
assumption had to be made that transitions from the 
ground state to higher excitations are already correctly 
described in the first few iterations. Accurate results 
in agreement with sum rules were obtained for the single 
particle spectrum in the absence of external fields. In the 
following, however, we demonstrate that this procedure 
is not rigorous and fails for the Anderson impurity model 
in a magnetic field. To remedy the defect, we introduce a 
new approach based on the concept of the reduced density 
matrix. This procedure (which in the following will be 
denoted as DM-NRG) makes use of the full information 
contained in iterative diagonalization and can therefore 



be considered as the true extension of Wilson's original 
work to dynamical quantities. 

To be specific, we consider the spin 1/2 Anderson 



Ho + iJimp where the impurity part is given 



model H 
by 



i^imp = V [JICq„ + /i.e. j + U Hf^Ufi - efUf-h S]. (1) 

Here we have introduced a local magnetic field h coupled 
to the impurity spin Sj, an on-site Coulomb repulsion U, 
and a hybridization A — ttV"^ /2 to the conduction band 
orbital cqo-. Units are chosen as Ti = ks = g — fJ-B — 1- 
Depending on the energy of the impurity level, e^-, differ- 
ent physical behaviour is realized. In the following, we 
focus on the symmetric (e/ = —U/2) and mixed valence 
(je^^l « A) regimes. The conduction band (extending in 
the range [—1,1]) is already written in the linear chain 
representation characteristic for NRG 
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(2) 



where the hopping matrix elements decay exponentially 
S,n ~ A~"/^ due to a logarithmic discretization of the 
conduction band. This model - while still a nontrivial 
many-body problem - can now be solved by iterative di- 
agonalization, keeping in each step only the lowest, most 
relevant levels. The number of iterations then corre- 
sponds to the temperature one is interested in accord- 
ing to T/v = cA~(^~^^/^, where c is a constant of order 
one. For calculating static quantities, all necessary infor- 
mation is thus obtained because only excitations on the 
scale T/v are relevant. For dynamical properties, however, 
an additional energy scale is provided by the frequency u 
which may be much larger than the temperature. Let us 
focus on the spin resolved single particle spectral density 



A, {UJ)=Y,\ HPM \^H^- E,n + ^n) 



(3) 



in the Lehmann representation where the |n) are the 
many-particle eigenstates of H and Z is the partition 



1 



function. Obviously, spectral information at frequen- 
cies Lo ^ T/v requires matrix elements between low- 
lying states and excitations which in iteration TV are 
not available anymore (they have already been lost by 
truncation). To circumvent this difficulty, the following 
"ad hoc" procedure has been used so far: In calculating 
A{u;), expression (|^) was simply evaluated in iteration 
step N' <^ N where T/v « w, assuming that for an anal- 
ysis of this spectral regime the low energy levels were 
described well enough. There is no rigorous argument 
to justify this assumption, as, for example, the crossover 
to the strong coupling fixed point and the corresponding 
change in the ground state may occur at a much lower 
temperature scale Tk ^ TV' . 
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FIG. 1. Comparison of single particle spectral functions 
for the symmetric model (A = 0.01, U = 0.1, e/ = -0.05) 
obtained by the method previously used ("NRG") and the 
generalized procedure presented here ("DM-NRG"). A small 
magnetic field h = 0.001 has been applied to the impurity. 

In fig. |l] we present results for the symmetric model (|^) 
at T = which have been calculated in this way. With- 
out external field, one obtains the well-known three-peak 
structure characteristic for a small Kondo temperature 
Tk. Switching on a small magnetic field h = 0{Tk) 
only affects the quasiparticle peak, while the high en- 
ergy spectrum is almost unchanged. This result is easily 
understood: In the iterations where the atomic levels 
are calculated, the NRG procedure does not yet "know" 
about the tiny perturbation that eventually breaks the 
spin symmetry of the ground state. One can, however, 
easily verify that this result is incorrect: Calculating the 
ground-state magnetization m (a static quantity) directly 
as a thermodynamic expectation value ((n/t ~^i/i)) and 
comparing with the value derived from the spectrum 

/o ^0 
dojA^iuj)- / duj Ai{uj) (4) 
-OC J — OQ 

the results do not agree (see table in fig. ^). Physically, 
the strong polarization of the impurity due to a small 
magnetic perturbation should suppress the upper atomic 
level because no particle excitations are possible any- 
more. This suppression is drastically underestimated by 



the method used so far (indeed, in the limit of vanishing 
Kondo temperature Tk it will not be seen at all). 

In order to capture this effect it is clearly necessary 
to obtain the correct ground state before calculating the 
spectra. This is achieved by the following two-stage pro- 
cedure: 

1) NRG iterations are performed down to the temper- 
ature T/v of interest, in particular we choose T/v <§; Tk to 
calculate ground-state properties. In each iteration step, 
we keep the information on the transformation between 
one set of eigenstates and the next, i.e. we save the cor- 
responding unitary matrix. After obtaining the relevant 
excitations at temperature T/v one can define the density 
matrix 

/5 = ^e"-^"/^" |m)/v(m| (5) 

rn 

which completely describes the physical state of the sys- 
tem. In particular, the equilibrium Green's function can 
be written as 

GT(t)=*0(i)Tr(p{/^(t),/|(O)}) (6) 
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FIG. 2. Reduced density matrix obtained by tracing out 
"environment" degrees of freedom of the chain. 

2) Now we repeat the iterative diagonalization for the 
same parameters. Each iteration step N' yields the 
single-particle excitations (and matrix elements of /^) 
relevant at a frequency lo ~ T/v/. But instead of using 
(^), we now employ (^ and evaluate the spectral func- 
tion with respect to the correct reduced density matrix 
[Q : As depicted in fig. ||, the complete chain is split into 
a smaller cluster of length N' and an environment con- 
taining the remaining degrees of freedom. In the product 
basis of these two subsystems, the full density matrix has 
the form 

/5 = ^ Pmini,m2n2l"^l)cnv|"l)sys('T'2KTO2| (7) 

which is in general not diagonal anymore. Performing a 
partial trace on the environment then yields the density 
submatrix 

E<tKWs(«2| (8) 
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(9) 



This projection is easily done using the previously stored 
unitary transformation matrices. Note that p'^°'^ - de- 
fined only on the shorter chain - contains all the rel- 
evant information about the quantum mechanical state 
of the full system. This concept has been applied very 
successfully in the density matrix renormalization group 
(DMRG) where the projection on a smaller subsys- 
tem is essential for obtaining eigenstates of the model. 
In NRG, on the other hand, diagonalization can be per- 
formed directly due to the logarithmic discretization, but 
to describe the efi^ect of the chain degrees of freedom on 
the impurity (or a small cluster) one has to determine 
p'^°'^. In the following, we therefore refer to the calcula- 
tional scheme presented here as DM-NRG. 

In fig. |l] we compare the spectrum calculated by the 
DM-NRG to the one obtained with the NRG version 
used so far in the literature (the same number of lev- 
els has been used in both calculations). The strong shift 
of spectral weight due to the polarized impurity is now 
clearly seen, as well as a slight change in the height and 
shape of the quasiparticle peak. The magnetization has 
been calculated from (^) for different values of h and is 
in good agreement with the static calculation (see fig. H). 
The remaining deviation of about three percent is due to 
an error in the total spectral weight. 
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FIG. 3. Impurity magnetization obtained by different 
methods: from the spectrum (ttinrg vs. mDM-NRc) and as 
a thermodynamic expectation value (mdircct). Impurity pa- 
rameters are chosen as A = 0.01 and U = 0.1. 




0.02 0.04 0.06 0.08 



FIG. 4. Shift of the spectral function with increasing mag- 
netic field at zero temperature. The impurity parameters are 
chosen as A = 0.01 and U = 0.1. 



The resulting field dependence of the spectrum is dis- 
played in fig. ^ With increasing /i, the Kondo reso- 
nance is suppressed and eventually merges with the lower 
atomic level. Regarding the total density of states (DOS) 
A{uj) = ^„ the Kondo peak is split by the field 

and the DOS at the Fermi level strongly reduced. This 
effect has been observed directly in measurements of the 
differential conductance through a quantum dot |2]. 



T=0.0055 
0.011 
0.022 
0.044 




FIG. 5. Temperature dependence of the spectrum at 
A = 0.1, U = 1.0 and h = 0.04. Note that at frequencies 
far below the temperature NRG does not yield any informa- 
tion. In this region the curves are fitted. 

So far calculations have been at T = 0. Upon increase 
of the temperature at a finite magnetic field, we expect 
a reduction of the average impurity magnetization due 
to thermal fluctuations. As a consequence, particle ex- 
citations with polarization in the field direction should 
gain spectral weight. In fig. ^ this effect is obvious: At 
temperatures T > /i, the asymmetry in v4|(lj) is strongly 
reduced. Note that in finite temperature NRG calcula- 
tions, no spectral information can be obtained at frequen- 
cies uj <S^T. In this region data have to be fitted. This 
important fact will be discussed in detail in a subsequent 
publication. 
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FIG. 6. Spin dependent spectral density at zero temper- 
ature for the asymmetric impurity with A = 0.01, U = 0.1 
and £/ = -0.02. 

Results for an asymmetric impurity close to the mixed 
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valence regime are shown in fig. |6[ The almost complete 
shift of spectral weight to the particle (hole) sector is 
again observed for the two spin polarizations, which in 
this case are not symmetric anymore. In the total DOS 
(fig. 1^), changes are less prominent. We merely observe a 
suppression of the quasiparticle peak and a redistribution 
of the corresponding weight to higher frequencies. 

Comparing our findings to previous calculations, it 
should be pointed out that up to now only the modi- 
fied perturbation theory ||T^ and the Quantum Monte 
Carlo method (QMC) have been applied to calculate 
the impurity spectrum in a magnetic field. The former is 
limited to small U, while QMC calculations have so far 
only been done in the mixed valence regime (and at tem- 
peratures T > Tk) due to the increase in computing cost 
for the symmetric case. In a recent NRG calculation on 
the Kondo model [|l8|, the problems discussed here did 
not occur due to the absence of atomic levels. Apart from 
these restrictions, we find qualitative agreement with our 
results, which do not suffer from similar limitations. 



broken spin symmetry. Our spectral results are in excel- 
lent agreement with the sum rule provided by the (static) 
magnetization. In the total density of states we find the 
splitting and suppression of the quasiparticle peak which 
is also observed experimentally. 

Future apphcations of the DM-NRG include DMFT 
calculations for phases with long range order, where 
symmetry-breaking perturbations and their effect on the 
spectrum have to be treated reliably. In addition, more 
complex impurity systems including orbital degeneracy 
may be studied, which (due to the rapid advances in 
nanoscale preparation techniques [Q) are of growing ex- 
perimental interest. 

The author would like to thank D. VoUhardt, R. Bulla 
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FIG. 7. Total spectral density Atot = Ai+Ai at zero tem- 
perature for the asymmetric impurity with A = 0.01, U = 0.1 
and 6/ = —0.02. Note that upon increasing h part of the spec- 
tral weight is shifted to the upper atomic level (not shown). 
The total weight is constant with high accuracy. 

In conclusion, we have presented a new method for cal- 
culating dynamical properties at arbitrary temperatures 
within the numerical renormalization group. It has been 
demonstrated that - despite logarithmic discretization 
- energy scale separation is in general not valid in the 
case of spectral quantities. This effect is neglected in the 
NRG scheme used so far in the literature. Within our 
generalized procedure (DM-NRG), based on the reduced 
density matrix, we can now account for changes in the 
ground state occuring at energies far below the external 
frequency scale. 

The method introduced here has been applied to the 
Anderson impurity in an external magnetic field, which 
is of great interest in view of recent transport measure- 
ments of quantum dots. Nonperturbative T = stud- 
ies have not been performed so far, mainly because of 
technical difficulties in extending NRG to systems with 
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